The Mott insulator LaTiOs in heterostructures with SrTiC>3 is metallic 
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It is shown that LaTiOs in superlattices with SrTiOa is not a Mott insulator but a strongly 
correlated metal. The tetragonal lattice geometry imposed by the SrTiC>3 substrate leads to an 
increase of the Ti 3d <2 S band width and a reversal of the t2 g crystal field relative to the orthorhombic 
bulk geometry. Using dynamical mean field theory based on finite-temperature multi-band exact 
diagonalization we show that, as a result of these effects, local Coulomb interactions are not strong 
enough to induce a Mott transition in tetragonal LaTi03. The metalicity of these heterostructures 
is therefore not an interface property but stems from all LaTi03 planes. 
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There exists currently considerable interest in the de- 
sign of nano-materials with electronic properties that dif- 
fer qualitatively from those of the constituent compo- 
nents in their bulk form. A particularly intriguing exam- 
ple is the formation of a thin metallic layer at the inter- 
face between two insulators. For instance, SrTiOs is a 
band insulator with an empty d band, whereas LaTiOs, 
with one d electron per site, is regarded as a textbook 
Mott insulator because of strong local Coulomb interac- 
tions Nevertheless, the pioneering work by Ohtomo 
et al. [2j shows that LaTi0 3 /SrTi0 3 superlattices are 
metallic, where the overall conductivity depends on the 
thickness of the LaTiC>3 interlayers and on the spac- 
ing between them. This phenomenon appears to follow 
from the fact that the Ti interface layer between adja- 
cent La and Sr planes formally exhibits a 3c? 5 valency, 
in contrast to the 3d 1 and 3d configurations of bulk 
LaTiOs and SrTiOs, respectively. Other examples are 
LaAlOs/SrTiOs heterostructures where a conducting in- 
terface was observed although both constituents are wide 
band gap perovskite insulators Q. Very recently, vari- 
ous other perovskite superlattices have been investigated 
experimentally @, S, @, 0, i ©, E3, EH • 

One of the hallmarks of transition metal oxides is their 
extreme sensitivity to small changes of key parameters 
such as temperature, pressure, impurity concentration, 
or structural distortion The aim of this work is to 
demonstrate that, as a result of a subtle change of crys- 
tal structure, LaTiOs in superlattices with SrTi03 is a 
strongly correlated metal, rather than a Mott insulator 
like bulk LaTiOs. As a consequence, the metalicity ob- 
served in Q is not only a property of the interface but of 
the entire LaTiOs layer. 

The origin of this symmetry induced insulator to 
metal transition is that, when LaTiOs is grown on cu- 
bic SrTiOs, the substrate imposes a tetragonal geometry 
on the first few layers EH • As we show below, this struc- 
tural modification leads to two significant changes of the 
electronic properties. First, it causes a substantial in- 
crease of the Ti tig band width. Second, the t 2g crystal 
field is weaker and has the opposite sign compared to bulk 



LaTiOs. Both effects ensure that the Mott transition in 
tetragonal LaTiOs occurs at a significantly larger criti- 
cal Coulomb energy U c than in bulk LaTiOs. Thus, for 
realistic values of U, thin LaTiOs layers in heterostruc- 
tures with SrTiOs are strongly correlated metals. These 
results underline the importance of carefully characteriz- 
ing the interfaces of superlattices that might be employed 
in future devices. 

The electronic properties of LaTiOs/SrTiOs het- 
erostructures have been studied theoretically within a 
variety of single-electron and many-electron models [HI, 
3, EE ES Eil, El 0, HI EH- The conditions under 
which LaTiOs in these superlattices might be a Mott in- 
sulator, however, have not yet been explored. In view of 
the striking interplay between Ti 3d ti g orbital degrees of 
freedom and local Coulomb interactions this issue is cru- 
cial and might differ considerably from the situation in 
bulk LaTiOs. As shown by Pavarini et al. the noncu- 
bic octahedral distortions in orthorhombic bulk LaTiOs 
give rise to nondiagonal components of the ti g density 
of states. Nevertheless, from these orbitals a new basis 
can be constructed (denoted here as a 9 , e' g ), in which the 
local density of states is nearly diagonal, and where the 
a g contribution lies about 200 meV below the nearly de- 
generate pair of e' g components, in qualitative agreement 
with experimental findings [23| . The local Coulomb en- 
ergy greatly enhances this Ti 3d orbital polarization so 
that for U « 5 eV a Mott transition occurs where the 
a g band is nearly half-filled and the e' g states are nearly 
empty [22|. The important conclusion from this picture 
is that the metal insulator transition in LaTiOs is inti- 
mately connected to the sign and magnitude of the ti g 
crystal field splitting which reflects, in turn, the orienta- 
tion and magnitude of the noncubic lattice distortion. 

Here we combine ab initio electronic structure calcu- 
lations for the SrTiOs induced tetragonal geometry of 
LaTiOs with finite-temperature dynamical mean field 
theory (DMFT) 24 , 2|| and compare the correlated elec- 
tronic properties with those of the bulk orthorhombic 
structure. As impurity solver we use a recent multi- 
band extension [2y| of exact diagonalization (ED) [27| 
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FIG. 1: (Color online) Upper panel: t^g density of states 
components LaTi03 in hypothetical tetragonal bulk geome- 
try. Solid (red) curve: d xy states, dashed (blue) curves: d xz , yz 
states. Ep = 0. Lower panel: analogous density of states 
components for central Ti atom in (LaTi03)6/(SrTiOs)4 su- 
perlattice (see text). 



which has been shown to be highly useful for the study 
of stron g c orrelations in several transition metal ox- 
ides [1^, HI, m, [3(|. Since this approach does not suf- 
fer from sign problems the full Hund exchange including 
spin-flip and pair-exchange contributions can be taken 
into account. Also, relatively large Coulomb energies and 
low temperatures can be reached. 

The structures of LaTi03 in the tetragonal bulk and 
superlattice geometries were optimized via the Vienna 
Ab-initio Simulation Package (VASP) [3l|, an imple- 
mentation of the projector augmented- wave method [321 ] . 
where the total energy of the system was calculated 
within non-spin-polarized GGA. The in-plane lattice con- 
stant was fixed to that of SrTiC-3, a = b = 3.92 A, 
whereas the lattice parameters in the z direction were 
fully optimized. We assumed a lxl geometry without 
considering rotation and tilting of TiOe octahedra. For 
each optimized geometry, we performed an LDA calcu- 
lation using a home-made FLAPW code to calculate the 
partial density of states projected on Ti 3d orbitals in 
muffin-tin spheres. For the bulk structure, the z axis 
layer spacing c = 4.01 A is only slightly smaller than 
c = 4.106 A obtained for a tetragonal lattice with the 
same unit cell volume as in the orthorhombic case. 




FIG. 2: (Color online) Density of states components of or- 
thorhombic bulk LaTiOa in nearly diagonal subband repre- 



sentation [221 ] 
green) curves: 



Solid (red) curve: a g 
e'„ states. Ef = 0. 



states, dashed (blue and 



Fig. 1 shows the Ti i 2 cj density of states components 
for the hypothetical tetragonal bulk geometry of LaTi03 
(upper panel). In this symmetry, the singly degenerate 
d xy density differs from the doubly degenerate d xz ^ yz den- 
sities. In the experimental work of Ohtomo et al. 0], 
rather thin LaTiC>3 layers containing up to five La planes 
were studied. The lower panel of Fig. 1 shows the t^g 
density of states components for a (LaTi03)e/(SrTi03)4 
superlattice, containing five Ti layers in a nominal d 1 
configuration (only La neighboring planes), three d° Ti 
layers (only Sr neighboring planes), and two d 5 Ti in- 
terface layers (La and Sr neighboring planes). Plotted 
is the density of states for the central Ti atom within 
the LaTiC>3 slab. The vertical height of the 06 octahe- 
dron surrounding this atom, 3.98 A, is slightly smaller 
than that in the tetragonal bulk, 4.01 A. As a result, the 
difference between the d xy and d xzyz densities is even 
smaller than that in bulk case shown in the upper panel 
of Fig. 1. 

The overall width of the tn g bands for both tetragonal 
structures is about Wtetra ~ 2.5 eV, where we have ne- 
glected small additional spectral weight in the region of 
the empty e g bands at higher energies. Also, the centroid 
of the singly degenerate d xy band is seen to lie above that 
of the d xz yz bands. As a result of this positive crystal 
field splitting, the subband occupancies (per spin band) 



are n T 



0.14 
0.15, n. 



n xz . yz ~ 0.18 for the tetragonal bulk and 
xz,yz ~ 0.175 for the tetragonal slab. 

*, g and e' g den- 



For comparison we show in Fig. 2 the a g 
sities of orthorhombic bulk LaTi03 [22| which are ob- 
tained from a linear transformation of the original Ti 
d xy , xz , yz orbitals. As mentioned above, in the a g , e' g ba- 
sis the orthorhombic density of states is nearly diagonal 
and the centroid of the a g density lies about 200 meV be- 
low the two e' g components, implying a negative crystal 
field splitting between the singly and doubly degener- 
ate densities. The subband occupancies (per spin band) 
are n a « 0.23 and n e / « 0.135. Thus, orbital po- 
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FIG. 3: (Color online) Ti 3d subband occupancies of LaTiO;i 
as a function of local Coulomb energy. Solid (red) curves: or- 
thorhombic bulk geometry, indicating metal insulator transi- 
tion at U » 5 eV. Short-dashed (blue) curves: tetragonal bulk 
geometry, long-dashed (green) curves: tetragonal slab geom- 
etry, both indicating Mott transitions at U > 6 eV. Since the 
crystal field splitting has opposite sign in the orthorhombic 
and tetragonal structures, the orbital polarization between 
singly degenerate a g (d xy ) and doubly degenerate e' g (d XZtVZ ) 
states is reversed. 



larization is much larger than for the tetragonal struc- 
tures shown in Fig. 1 and has the opposite sign. More- 
over, the width of the a g , e' g bands, W or tho ~ 2 eV, is 
considerably smaller than in the tetragonal case. On 
the whole, therefore, the noncubic distortions are sig- 
nificantly smaller for the tetragonal interlayer structure 
than for the orthorhombic bulk form of LaTi03. This 
difference is also apparent in the shape of the density of 
states since the tetragonal distributions shown in Fig. 1 
resemble more closely the characteristic ti g density of cu- 
bic SrV0 3 (width W cuhic = 2.8 eV (H). 

We now show that the combined influence of these dif- 
ferences: the wider band width, the weaker orbital polar- 
ization, and the opposite sign of this polarization, implies 
that LaTiOs in the tetragonal superlattice geometry is a 
strongly correlated metal, rather than a Mott insulator 
as in its orthorhombic bulk form. 

To evaluate the electronic properties of LaTiOs in the 
presence of strong local Coulomb interactions we use 
multi-band ED/DMFT [26]. Full Hund exchange is in- 
cluded, with the exchange energy J = 0.65 eV held fixed 
at the value appropriate for the bulk material j2^|. To 
locate the Mott transition, the intra-orbital Coulomb en- 
ergy U is varied, and the inter-orbital Coulomb inter- 
action is given by U' = U — 2 J. The temperature is 
T = 20 meV. 
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FIG. 4: (Color online) LaTiOs 3d t2 g spectral distributions. 
Upper panel: tetragonal superlattice geometry, U = 5.5 eV, 
indicating strongly correlated metallic behavior. Lower panel: 
orthorhombic bulk geometry, {7 = 5 eV, indicating Mott in- 
sulating behavior. Solid (red) curves: singly degenerate d xy 
(a g ) states, dashed (blue) curves: doubly degenerate d xz>yz 
(e' g ) states; Ef = 0. 



suits using the Quantum Monte Carlo technique as im- 
purity solver [22j. The spectral distributions (see Fig. 4) 
show that the nearly half-filled a g band then consists of 
lower and upper Hubbard peaks whereas the e' g bands are 
pushed above the Fermi level 22], 3(| ■ In the case of the 
tetragonal structures, the situation is qualitatively dif- 
ferent since the crystal field is weaker and has opposite 
sign. Now the enhancement of the orbital polarization 
via the Coulomb interaction implies that the d xy band 



xy 

is pushed above the Fermi level and that the degenerate 
dxz.yz bands are driven to 1/4 occupancy. As a result 
of this reversal of crystal field, and because of its weaker 
magnitude (i.e., weaker orbital polarization in the un- 
correlated limit), the metal insulator transition occurs 
above U = 6 eV. Since the local Coulomb energy for Ti 
is approximately 5 eV [22j , this result implies that in the 
tetragonal geometry LaTiOs is on the metallic side of 
the Mott transition. As mentioned above, we have ne- 
glected small contributions to the tetragonal density of 
Fig. 3 shows the a g , e' subband occupancies as a func- states in the re S ion of the e a bands - Inclusion of 
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tion of Coulomb energy. For the orthorhombic geom- 
etry, nearly complete orbital polarization is reached at 
about U = 5 eV, which agrees with previous DMFT re- 



these states would increase the tetragonal ti g band width 
and reinforce the metallic behavior. 

As can be seen in Fig. 3, the Mott transition for the 



4 



tetragonal slab geometry occurs at slightly larger U than 
for the tetragonal bulk. This shift is a consequence of 
the smaller orbital polarization in the former case, asso- 
ciated with the weaker deviations of the density of states 
components from cubic symmetry, as is evident also from 
the results shown Fig. 1. 

In Fig. 4 we compare the LaTiC>3 t^g spectral distri- 
butions for the tetragonal slab and orthorhombic bulk 
crystal structures. Since we are concerned here only with 
the difference between metallic and insulating behavior, 
it is sufficient to use the ED cluster spectra which can be 
evaluated at real u>, without requiring extrapolation from 
imaginary Matsubara frequencies. In the orthorhombic 
case, the system is already insulating at U = 5 eV, with 
an excitation gap formed between the lower Hubbard 
peak of the a g band and the empty e' g bands, in agree- 



ment with Ref. [221 ■ The tetragonal slab geometry, on 
the other hand, is still metallic at U = 5.5 eV, with con- 
duction states near Ep from d xy and d xz ^ yz bands. 

In the DMFT calculations discussed above we have 
focused on the central Ti plane in the tetragonal LaTiC>3 
layer. Ti planes at the interface with SrTiC>3 are also 
metallic because of the nominal 3d 5 occupancy 
The mutual enhancement of both effects ensures that the 
entire LaTiC>3 layer is metallic. 

The above results demonstrate the extreme sensitiv- 
ity of the electronic properties of perovskites to small 
changes in key parameters. In the present case, the 
slightly modified crystal structure caused by the depo- 
sition of LaTiC>3 onto SrTiC>3 gives rise to a subtle, but 
important change in its single-electron spectral distribu- 
tion, which, in turn, has a profound effect on the prop- 
erties of the correlated electron system. For an adequate 
interpretation of experimental data it is therefore of great 
importance to have accurate structural information. On 
the theoretical side, it is crucial to include the full details 
of the single-particle electronic properties for a given lat- 
tice geometry and to account for the orbital degrees of 
freedom in the multi-band many-body calculation. 

In summary, we have shown that in LaTiOs/SrTiOs 
superlattices the change from orthorhombic to tetrago- 
nal crystal structure enforced by the SrTiC>3 substrate 
causes a fundamental change in the electronic proper- 
ties of LaTiC>3 . Instead of being a Mott insulator as in 
the bulk geometry, it is now a highly correlated metal. 
Thus, the metalicity observed in these heterostructures 
is not only associated with the interfaces but with the 
entire LaTiC>3 layers. It would be very interesting to 
use transmission electron microscopy to determine the 
LaTiC>3 layer thickness at which the orthorhombic struc- 
ture begins to be energetically favorable. Mott insulat- 
ing behavior should then emerge in the central regions of 
these LaTiC>3 layers at sufficiently large thickness. 
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